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Abstract 

We present a systematic method for computing explicit approximations to martingale 
representations for a large class of Brownian functionals. The approximations are obtained 
by obtained by computing a directional derivative of the weak Euler scheme and yield 
a consistent estimator for the integrand in the martingale representation formula for any 
square-integrable functional of the solution of an SDE with path-dependent coefficients. 
Explicit convergence rates are derived for functionals which are Lipschitz-continuous in the 
supremum norm. Our results require neither the Markov property, nor any differentiabil¬ 
ity conditions on the functional or the coefficients of the stochastic differential equations 
involved. 
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1 Introduction 


Let W be a standard d-dimensional Brownian motion defined on a probability space (fl, J 7 , IP) and 
(Ft) its (P-completed) natural filtration. Then, for any square-integrable J-r-measurable random 
variable H , or equivalently, any square-integrable (J-^-martingale Y(t) = K[H\Ft], there exists 

a unique (J 7 t )-predictable process <f) with E / tr(< ^(u) 1 <j>(u))du 


< oo such that: 


H = Y(T) = E [H] + 


<t> ■ dW. 


( 1 ) 


The classical proof of this representation result (see e.g. [33]) is non-constructive. However 
in many applications, such as stochastic control or mathematical finance, one is interested in 
an explicit expression for </>, which represents an optimal control or a hedging strategy. Ex¬ 
pressions for the integrand </> have been derived using a variety of methods and assumptions, 
using Markovian techniques JO. |T3] US] [23 , 30] , integration by parts [2] or Malliavin calculus 
E E3 M M M Hi ESI- Some of these methods are limited to the case where Y is a Markov 
process; others require differentiability of H in the Frechet or Malliavin sense n ni ns na or 
an explicit form for the density [2]. Almost all of these methods invariably involve an approx¬ 
imation step, either through the solution of an auxiliary partial differential equation (PDE) or 
the simulation of an auxiliary stochastic differential equation. 


A systematic approach to obtaining martingale representation formulae has been proposed 
in jj], using the Functional Ito calculus B21E1EI: it was shown in [7] Theorem 5.9] that for any 
square-integrable (J 7 t)-martingale Y, 


Vf<E[0,T], Y(t) = Y(0)+ [ V w Y-dW P-a.s. 

Jo 

where is the weak vertical derivative of Y with respect to W, constructed as an L 2 limit 

of pathwise directional derivatives. This approach does not rely on any Markov property nor on 
the Gaussian structure of the Wiener space and is applicable to functionals of a large class of Ito 
processes. 


In the present work we build on this approach to propose a general framework for computing 
explicit approximations to the integrand <j> in a general setting in which H is allowed to be 
a functional of the solution of a stochastic differential equation (SDE) with path-dependent 
coefficients: 

dX(t) = b(t, X t )dt + cr(t, X t )dW{t) X(0) = x 0 G R d (2) 

where X t = X(t A .) designates the path stopped at t and 

b : [0, T] x L>([0,T],R d ) ->• R d , a: [0,T] x D([0,T],R d ) ->• M d (R) 

are continuous non-anticipative functionals. For any square-integrable variable of the form H = 
g(X(t), 0 <t<T ) where g : (D([0, T], K d ), ||.||oo) — > M is a continuous functional, we construct 
an explicit sequence of approximations <j> n for the integrand <f> in 0 . These approximations 
are constructed as vertical derivatives, in the sense of the functional Ito calculus[7], of the weak 
Euler approximation F n of the martingale Y, obtained by replacing X by the corresponding 
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Euler scheme n X : 


(t>n{t) = X bJ F n (t,W), where F„(t,w) = E 


5(n^(w © WO) 


where © is the concatenation of paths at t and V^F n is the Dupire derivative 01 a direc¬ 
tional derivative defined as a pathwise limit of finite-diference approximations, and thus readily 
computable path-by-path in a simulation setting. 


The main results of the paper are the following. We first show the existence and continuity 
of these pathwise derivatives in Theorem 14.11 The convergence of the approximations <p n to the 
integrand </> in © is shown in Proposition 15.11 Under a Lipschitz assumption on g , we provide 
in Theorem an L p error estimate for the approximation error. The proposed approximations are 
easy to compute and readily integrated in commonly used numerical schemes for SDEs. 


Our approach requires neither the Markov property of the underlying processes nor the differ¬ 
entiability of coefficients, and is thus applicable to functionals of a large class of semimartingales. 
By contrast to methods based on Malliavin calculus m is m m m [18], it does not require 
Malliavin differentiability of the terminal variable H nor does it involve any choice of ’Malliavin 
weights’, a delicate step in these methods. 


Ideas based on Functional Ito calculus have also been recently used by Leao and Ohashi 
[25] for weak approximation of Wiener functionals, using a space-filtration discretization scheme. 
However, unlike the approach proposed in ESI, our approach is based on a Euler approximation 
on a fixed time grid, rather than the random time grid used in [26], which involves a sequence of 
first passage times. Our approach is thus much easier to implement and analyze and is readily 
integrated in commonly used numerical schemes for approximations of SDEs, which are typically 
based on fixed time grids. 


Outline We first recall some key concepts and results from the Functional Ito calculus in section 
[2] Section [3] provides some estimates for the path-dependent SDE © and studies some properties 
of the Euler approximation for this SDE. In SectionUwe show that the weak Euler approximation 
(Definition [9]) may be used to approximate any square-integrable martingale adapted to the 
filtration of X by a sequence of smooth functionals of X, in the sense of the functional Ito calculus. 
Moreover, we provide explicit expressions for the functional derivatives of these approximations. 
Section [5] analyzes the convergence of this approximation and provides error estimates in Theorem 
15.11 Finally, in Section Owe compare our approximation method with those based on Malliavin 
calculus. 

Notations: In the sequel, we shall denote by -A/l f /. rl (M) the set of all d x n matrices with real 
coefficients. We simply denote R d = At^i^R) and .Ad^R) = A4d,d(R)- For A £ A4rf(M), we shall 
denote by t A the transpose of A, and ||H|| = \J tr (^AA) the Frobenius norm of A. For i,i/g R d , 
x ■ y is the scalar product on R d . 

Let T > 0. We denote by D([0, T], R d ) the space of functions defined on [0, T] with values in R d 
which are right continuous with left limits (cadlag). For a path lo £ D([0,T],M d ) and t £ [0, T], 
we denote by: 

• u(t) the value of ui at time t, 
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• uj(t —) = lim w(s) its left limit at t , 

s—>t,s<t 

• tut = ai(t A •) the path of w stopped at t 

• w t- = W 1[CM) + u(t— )l[t,T] 

• ||w||oo = sup{|w(f)|, f £ [0,X 1 ]} the supremum norm. 

We note that w t and co t - are elements of Z?([0, T], R d ). For a cadlag stochastic process X , we 
shall similarly denote X t (.) = X(t A .) and X t - = Xl[ 0jt ) + X(t— )l[t,T]- 


2 Functional Ito calculus 


The Functional Ito calculus is a functional calculus which extends the Ito calculus to path- 
dependent functionals of stochastic processes. It was first introduced in a pathwise setting 
nanaiia using a notion of pathwise derivative for functionals on the space of right-continuous 
functions with left limits, and extended in [7] to a weak calculus applicable to all square-integrable 
martingales, which has a natural connection to the martingale representation theorem. We recall 
here some key concepts and results of this approach, following [1] . 

Let X be the canonical process on SI = Z?([0,T],R d ), and be the filtration 

generated by X. We are interested in non-anticipative functionals of X , that is, functionals 
F : [0,T] x Z?([0, T],R d ) H> R such that 

Vw S fl, F(t,cu) = (3) 

The process t >->• F(t,X t ) then only depends on the path of X up to t and is (.F t 0 )-adapted. 

It is convenient to define such functionals on the space of stopped paths [4]: a stopped path 
is an equivalence class in [0,T] x D([0,T],R d ) for the following equivalence relation: 

(t,uj) ~ (t',u/) (i = t' and uj t =uj , t ,). (4) 

The space of stopped paths is defined as the quotient of [0, T] x U([0, T], R d ) by the equivalence 
relation 

A t = {(t, ui(t A •)), (t,u) € [0, T] x £([0,T],M d )} = ([0,T] x D([0,T],R d )) / ~ 

We denote Wt the subset of At consisting of continuous stopped paths. We endow this set with 
a metric space structure by defining the following distance: 

doo((t,u>), (t',u>')) = sup \uj(u A t) - ui'{u A t')\ + \t - t'\ = ||w t - u' t , ||oo + \t - t'\ 
ue[o,T] 

(At, doo) is then a complete metric space. Any functional verifying the non-anticipativeness 
condition m can be equivalently viewed as a functional on F : At —» R: 

Definition 1. A non-anticipative functional on D([0, T], R d ) is a measurable map 
F: (At, doc) —> R on the space (At,cIoo) of stopped paths. 
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Using the metric structure of (Aj^doo), we denote by C 0,0 (A'r) the set of continuous maps 
F : (ATjdoo) i —y R. Some weaker notions of continuity for non-anticipative functionals turn out 
to be useful |J>]: 

Definition 2. A non-anticipative functional F is said to be: 


• continuous at fixed times if for any t £ [0,T], F{t,-) is continuous with respect to the 
uniform norm || • ||oo in [0,T], i.e. Vw £ D([0, T], R d ), Ve > 0, 3r/ > 0, Vu/ £ D([0, T], R d ), 

IK - u' t \\ < rj => \F(t,u) - F{t,u')\ < e 

• left-continuous ifV(t,u>) £ At, Ve > 0, > 0 such that M(t' ,uj') £ At, 

( t'<t and doo ((t, to), (t 1 , to 1 )) < rj) ==> \F(t, w) — F(t', u/)| < e 

We denote by C^’°(At) the set of left-continuous functionals. Similarly, we can define the 
set C®’°(At) of right-continuous functionals. 


We also introduce a notion of local boundedness for functionals. 

Definition 3. A non-anticipative functional F is said to be boundedness-preserving if for every 
compact subset K ofR d , Vfo £ [0, X 1 ], 3C(K,to) > 0 such that: 

Vte[0, to], V(t,ui) £ A t , w([0, t]) C K =>• F(t, w) < C{K, t 0 ). 

We denote by B(At) the set of boundedness-preserving functionals. 


We now recall some notions of differentiability for functionals following mm. For e G R d and 
w £ D([0, T], R d ), we define the vertical perturbation of (t, w) as the cadlag path obtained by 
shifting the path by e after t: 

lOt = uj t + el[ ti 7’]. 

Definition 4. A non-anticipative functional F is said to be: 


• horizontally differentiable at (t, u>) £ At if 


VF(t, ui) 


lim 
h —►()+ 


F(t + h,w) — F(t,oj ) 
h 


exists. If VF{t,ui) exists for all ( t,w) £ At, then the non-anticipative functional VF is 
called the horizontal derivative of F. 

• vertically differentiable at (t, ui) £ At if the map: 

R d —*• R 

e H > F(t,u>t + e\yt,T]) 

is differentiable at 0. Its gradient at 0 is called the vertical derivative of F at (t,u>): 
VuF(t,u) = ( diF(t,w),i = !,••■ , d) £ R d 
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with 


diF{t,uj) = lim 
h —>-0 


F(t,LJt + heil [t ^) - F(t,u)t) 


where (e*, i = 1, • • 
{t, w) £ A r , V W F 
of F. 


,d) is the canonical basis of R d . If F is vertically differentiable at all 
(■ t , u>) defines a non-anticipative map called the vertical derivative 


We may repeat the same operation on V W F and define similarly V^F, V^F, • • •. This leads 
us to define the the following classes of smooth functionals: 

Definition 5 (Smooth functionals). We define Ck’ fc (A-r) as the set of non-anticipative function¬ 
als F : (A t, d) —> R which are 


• horizontally differentiable with VF continuous at fixed times; 

• k times vertically differentiable with V^F £ C°'°(At) for j = 0, • • • ,k; 

• VF, VuF, • • ■ ,V*F£B(A r ). 

We denote Cj ,00 (At) = rifc>iCj’ fc (AT). 

Many examples of functionals may fail to be globally smooth, but their derivatives may still 
be well behaved except at certain stopping times, which motivates the following definition 0J: 

Definition 6. A non-anticipative functional F is said to be locally regular of class Cj’ (Ar) if 
there exists an increasing sequence (t„)„> o of stopping times with tq = 0 and r n —► oo, and a 

n—f oo 

sequence of functionals F n £ C^’ 2 (A t) such that: 

- F «(^ W ) 1 [TnM,r„ +1 (a;))(f), V(t,w) £ A t 

n> 0 


We recall now the functional fto formula for non-anticipative functionals of a continuous 
semimartingale [7, Theorem 4.1]: 

Proposition 2.1 ([B][7]). Let S be a continuous semimartingale defined on a probability space 
(ff,F,P). For any non-anticipative functional F £ C 1c ] c (At) and any t £ [0,T], we have: 

F(t, S t ) — F(0, So) = [ VF(u,S u )du+ [ V u F(u, S u ) ■ dS(u) + \ [ tr (V 2 F(u, S u )d[S](u)) 

Jo Jo z Jo 

Actually the same functional fto formula may also be obtained for functionals whose vertical 
derivatives are right-continuous rather than left-continuous. We denote by C^’ (At) the set of 
non-anticipative functionals F satisfying: 


• F is horizontally differentiable with VF continuous at fixed times; 

• F is twice vertically differentiable with F £ C°’°(A t) and V U F, V^F £ C°’°(At’); 
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. DF,V w F,V 2 Fg B(A r ); 


The localization is more delicate in this case, and we are not able to state a local version of the 
functional Ito formula by simply replacing F n G C^’ 2 (At) by F n G C^’ 2 (At) in Definition [6] 
(see Remark 4.2 in HU). However if the stopping times r„ are deterministic, then the functional 
Ito formula is still valid (Proposition 2.4 and Remark 4.2 in IB11. 

Definition 7. A non-anticipative functional is said to be locally regular of class Cj^^At) if 
there exists an increasing sequence (t n ) n > o of deterministic times with to = 0 and t n — > oo, 

n—too 

and a sequence of functionals F n G C^’ 2 (At) such that: 

F(t,u) = ^F„(t,u;)l [Wn+l) (t), V(t,w) G A t 

n> 0 

Proposition 2.2 ([7]). Let S be a continuous semimartingale defined on a probability space 
(Cl, .F,P). For any non-anticipative functional F G Cjo CJ .(At) and any t G [0,T], we have: 

F(t,S t )-F(0,S o ) = [ VF(u,S u )du+ [ V„F(u,S u ) ■ dS(u) 

Jo Jo 

tr (V 2 F(u, S»)d[S\ (u)) P-a.s. 

Finally we present briefly the martingale representation formula established in [7]. Let 
(-^t)te[o,T] be a continuous Revalued martingale defined on a probability space (Cl,F, P) with 
absolutely continuous quadratic variation: 

[X](t) = f A(u)du 

Jo 

where A is a A4d(R)-valued process. Denote by (F*) the natural filtration of X and Cl' 2 (X) 
the set of (Jy X )-adapted processes Y which admit a functional representation in C^’ (At): 

Cl’ 2 (X) = {Y, 3F G C^’ 2 (At), Y(t) = F(t, X t ) dt x dP-a.e.} (5) 

If Aft) is non-singular almost everywhere, i.e. det(A(f)) ^ 0, dt x dP-a.e., then for any Y G 
Cl’ 2 (X), the predictable process 


X x Y(t)=V ul F(t,X t ) 

is uniquely defined up to an evanescent set, independently of the choice of F G C^’ 2 (At) in the 
representation ©. This process VxY is called the vertical derivative of Y with respect to X. 
For martingales which are smooth functionals of X, the operator V.y : Cl' 2 (X ) i —> C ( 0,0 (X) yields 
the integrand in the martingale representation theorem: 

Corollary 2.1. IfYG Cl' 2 (X) is a square-integrable martingale, then 

Vf G [0, T\, Y(t) = Y(0)+ f VxY-dX P-a.s. 

Jo 
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Consider now the case where X is a square-integrable martingale. Let A4 2 (X) be the space 
of square-integrable ( F /)-martingales with initial value zero, equipped with the norm ||U|| 2 = 
Y / E|V(T)| 2 . Cont & Fournie j?, Theorem 5.8] show that the operator Vx : Cl’ 2 (X) i-»- C ; ° : °(A') 
admits a unique continuous extension to a weak derivative V.v : M 2 (X) — ► C 2 (X) which satisfies 
the following martingale representation formula: 

Proposition 2.3 ([7]). For any square-integrable (F/')-martingale Y, we have: 

VtG[0,T], Y{t) = Y{ 0)+ / XxY-dX P-a.s. 

do 

This weak vertical derivative XxY coincides with the pathwise vertical derivative X u F(t, X t ) 
when Y admits a locally regular functional representation, i.e. Y(t) = F(t,Xt) with F G 
U C^ 2 r (A r ). For a general square-integrable martingale Y, the weak derivative XxY 
is not directly computable through a pathwise perturbation. An approximation procedure is thus 
necessary for computing X \Y. The result of [7] guarantees the existence of such approximations; 
in the sequel we propose explicit, and computable, constructions of such approximations. 


3 Euler approximations for path-dependent SDEs 


Let IT be a standard d-dimensional Brownian motion defined on a probability space (Id, F, P) 
and (Ft) its (P-completed) natural filtration. We consider the following stochastic differential 
equation with path-dependent coefficients ©: 

dX(t) = b(t, X t )dt + aft, X t )dW(t), X(0) = x 0 G 

where b : A t —» R d , a : A t —> M^(K) are non-anticipative maps, assumed to be Lipschitz 
continuous with respect to the following distance d defined on A^: 

d((t,Uj),(t',U]')) = SUp \uj(u A t) - Uj'(u A t')\ + y/\t - t'\ = || iO t - Wj/||oo + V\ t ~ t'\ 

«e[o,T] 

Assumption 1 (Lipschitz continuity of coefficients), b : (A r,d) —> K d ,u : (A T,d) —> Mj(R) are 
Lipschitz continuous: 

3K Lip > 0, Vt,t' G [0,T],Vu;V G D([0,T},R d ), 

| b(t,uj) - b(t',J )| + ||er(t,w) - <j(t',u')\\ < K Lip d((t,<jj), (t\o/)). 

Remark 3.1. This Lipschitz condition with respect to the distance d is weaker than a Lipschitz 
condition with respect to the distance doo introduced in the previous section: it allows for a Holder 
smoothness of degree 1/2 in the t variable. 


Under Assumption Q] (J2]) has a unique strong, (J 7 t )-adapted, solution X. 

Proposition 3.1. Under Assumption |7J there exists a unique (F/j-adapted process X satisfying 
©. Moreover for p > 1, we have: 

E[||AVia <C(l + \x 0 \ 2 ne CT (6) 

for some constant C = C(p, T, Kn p ) depending on p , T and Kn p . 
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Remark 3.2. Assumption^ might seem to be quite strong. Indeed, the previous proposition still 
holds under weaker conditions. For example, the Holder condition with respect to differences in 
t can be replaced by the weaker condition: sup \b(t, 0| + ||cr(t, 0)|| < oo where 0 denotes the path 

te[o,T] 

which takes constant value 0. However, this assumption is necessary for the convergence of the 
Euler approximation described later in this section, and especially the results concerning its rate 
of convergence. 


Proof. Existence and uniqueness of a strong solution follows from 22] (Theorem 7 , Chapter 5): 
see {3J Section 5]. Let us prove ©. Using the Burkholder-Davis-Gundy inequality and Holder’s 
inequality, we have: 


E[l|Xr||£] < C(p) |* 0 | 2p + E 


< C(p,T) |x 0 | 2p + E 


\b(t, X t )\ 2 dt 


•E 


\<r(t,X t )\\ 2 dt 



' r T „ 

+ E 

/ \\<j(t,X t )\\ p dt 


Jo 


\b(t,X t )\^dt 

_ J 1 ^ 

< C{p,T)(\x 0 | 2p + E jI (16(0,0)1 + 11^(0,0)|| +K Lip {Vi+\\X t \\ 00 )y P dt 

< C(p,T,K Lip ) (\ Xo \ 2p + l + ^ E||X t ||£dtj 


And we conclude by Gronwall’s inequality. 


□ 


In the following, we always assume that Assumption [Tj holds. The strong solution X of 
equation ((2J is then a semimartingale and defines a non-anticipative functional X : Wt —> 
given by the Ito map associated to ©• 


3.1 Euler approximations as non-anticipative functionals 

We now consider an Euler approximation for the SDE © and study its properties as a non- 

T 

anticipative functional. Let n G N, 6 = —. The Euler approximation n X of X on the grid 

n 

(tj = jS,j = 0 ..n) is defined as follows: 

Definition 8. [Euler scheme] For oj G _D([0, T],R d ), we denote by n X(u) G D([0,T],K d ) the 
piecewise constant Euler approximation for © computed along the path to, defined as follows: 
n X(ui) is constant in each interval [tj,tj+ 1 ), V0<j < n — 1 with n X(0 ,ui) = Xo and 

n X(t j+1 ,uj) = n X(tj,u>) + b(tj,„X tj (u>))5 + a(tj, n X t .(ui))(u>(tj + 1~) - uj(tj-)), (7) 

where n X t (u>) = n X(t A ,,ui) and by convention w(0—) = w(0). 


When computed along the path of the Brownian motion W, n X(W) is simply the piecewise 
constant Euler-Maruyama scheme m for the stochastic differential equation ©• 
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By definition, the path n X{ui) depends only on a finite number of increments of u: u>(t\—) — 
w(0), • • •, co(t n —) — oj{t n - 1 —). We can thus define a map 

P«:-M d , n (M)^Z)([0,T],M d ) 

such that for w G D([0, T]), M. d ) 

Pn(v(ti-) - U)(0),U){t 2 ~) - - ■ ■ ,u(t n -) - u(t n - 1~)) = n X(w). (8) 

By a slight abuse of notation, we denote pt(y) the path p n (y) stopped at t. 

The map p n : Ald,n(®0 —► (D([0, T], R d ), || • ||oo) is then locally Lipschitz continuous, as 
shown by the following lemma. 

Lemma 3.1. For every p > 0, there exists a constant C(p, Kn p ,T) such that for any y = 
(yi,--- ,yn),y' = (y[,--- ,y' n ) e 7W d , n (K), 

max \y k \ V \y' k \ < p => \\p n (y) ~ Pn{y')\\oo < C(p, K Lip , T) max \y k - y' k \. 

1 <k<n 1 <k<n 


Proof. As the two paths p n (y) and p n (y') are stepwise constant by construction, it suffices to 
prove the inequality at times {tj)o<j< n - We prove by induction that: 

\\Ptj(y) -Pt,(y')lloo < C(p,K Lip ,T) max \y k - y' k \ (9) 

1<K<J 

with some constant C which depends only on p , I\n p ,T (and n). 

For j = 0, this is clearly the case as p(y)( 0) = p(y’)( 0) = Xq. Assume that © is verified for 
some 0 < j < n- 1, consider now || p tj+1 (y) ~ Pt j+1 {y') ||oo, we have: 

Pn{y){tj+ 1 ) =Pn{y){tj) +b(t j ,p tj (y))6 + a(tj,p tj (y))y j+1 


and 

Thus 


Pn(y') (tj+ 1 ) = Pn ( y ') (tj) + b(tj , Ptj ( y'))S + a(tj , Pt . (y'))y' j+1 ■ 


\Pn{y)(tj+l) - Pn(y')(tj+l)\ 

< \Pn(y)(tj) -Pn(y')(tj)\ + \b{tj,p tj {y)) - b{t h p tj {y'))\S 

+\W( t j,Pt j (y))\\ ■ \yj+i - y'j+i I + \W(tj,Ph (y)) - v(tj,Ptj(y'))\\ ■ \y' j+ il 

< C(p , K Lip , T ) max \y k - y' k \ + K Lip C{p , K Lip , T) max \y k - y' k \S 

l<k<j 1<«<7 

+ (||cr(0, 0)|| + K L i P {y/tj + |bt,.(y))||oo)) \yj +1 - y'j+i\ + K Lip C{p, K Lip ,T)p max |y fc - y' k 

< C(??, A'iip, T) max \y k - y' k \ 


(The constant C may differ from one line to another). 

And consequently we have: 

\\pt j+ 1 {y) -Pt j+1 {y')\\oo < C(p, K Lip , T) max \y k - y' k \ 

for some different constant C depending only on p, K k i P and T (and n). And we conclude by 
induction. □ 
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3.2 Strong convergence 


To simplify the notations, n X- t(Wt) will be noted simply n X T in the following. The following 
result, which gives a uniform estimate of the discretization error, Xt — extends similar 
results known in the Markovian case 14, 31, 22] to the path-dependent SDE (J2J (see also 120]): 

Proposition 3.2. Under Assumption]]^ we have the following estimate in L 2p for the strong 
error of the piecewise constant Euler-Maruyama scheme: 

E\ sup ||X(s)- n X(s)|| 2 M <C(x 0 ,p,T,K Lip ) f 1 + ^°S n Y V p > 1 

\«e[o,T] J V n ) 

with C a constant depending only on Xq, p, T and KLi p . 


Proof. The idea is to construct a ’Brownian interpolation’ n X T of n X T : 

n X(s)=xo+ b (u, n Xj du+ a (u, n XdW(u) 
Jo Jo 


where u = 


L J 


• 8 is the largest subdivision point which is smaller or equal to u. 


Clearly n X is a continuous semimartingale and || sup |X(s) — n X(s)||| 2 P can be controlled 

se[o,T] 

by the sum of the two following terms: 

|| sup |X(s) - „X(s)||| 2 p < || sup |X(s) - „X(s)||| 2 p + || sup |„X(s) - n X(s)||| 2 p (10) 

se[0,T] «6[0,T] s£[0,T] 

We start with the term || sup |X(s) —„W(s)||| 2 P . Using the Burkholder-Davis-Gundy inequality 

se[o,T] 

and Holder’s inequality, we have 


E||X t - n X T \\ 2p < C{p) E 


( 

’ r T 

2 P 

' j'T 2 

E 

/ I b(s,X s ) - b(s, n X s )\ds 

+ E 

/ cr(s, X s ) - (r(s,A) ds 

V 

Jo 


Jo 


< C(p,T)\E 


r T 

/ ||cr(s, X s ) — cr(s, n Xs) 11 ds 

Jo 


< C(p,T,K Lip )E 


\b(s,X s )-b(s, n X ± )\ 2p ds 

[ T ((s-s) p + ||X 8 - ds 

Jo 

< C(p,T, K L ip) (^ + j\\\X s - n X s \\%ds'\ 

We have used n X s = n X s as n X is piecewise constant. 

Consider now the second term || sup | n X(s) — n X(s)||| 2 P - Noting that: 

s6[0,T] 

n X{s) - n X{s) = n X{s) - n X{s) = b (a, (a - a) + a (a, (W(s) - W(s)), 
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we have 


||n-^T — nXj' II oo < C(KLi pi T)( 1 + Hji^tIIoo) 



sup |W(s) — W(s)| 

se[o,T] 


and 


iX T - n X T \\% < C(p,K Lip ,T)^ E[(l + || n X r || 0O ) 2 *’] 


+C(p,K Lip ,T) E 


(1 + IlnXrlloo) sup |W(s)-W(s)| 
«e[o,T] 


2 P 


By the Cauchy-Schwarz inequality, we have: 

n\nX T - n X T fZ < C(p, K L i p , T) (l + yj E|| n X T ||^) 
We will make use of the following result .29] p. 203]: 


E sup |W(s) — W'(s)| 4p 
«e[o,T] i 


( 11 ) 


Vp > 0, || sup \W(s)-W(s)\\\ p <C(W,p)J—(l+\ogri) 

s€[0,T] V n 

which results from the following lemma |29[ Lemma 7.1]: 

Lemma 3.2. Let Y\, ■ ■ ■ , Y n be non-negative random variables with the same distribution sat¬ 
isfying E (e AYl ) < oo for some A > 0. Then we have: 

Vp > 0, ||max(Yi, • • • ,F„)|| p < ^(logn + C(p, Yi, A)) 


We have thus: 


f 


sup 

«e[o,T] 


\W{s)-W{s)\ i P <C{p,T) 


1 + log n 


n 


( 12 ) 


Furthermore, using again the Burkholder-Davis-Gundy inequality, we have: 


E||»* r ||£ 


< E||A||£ 

< C(p)lxt P +E 


J \HS;nXs)\d.s\ 



< C(xo,p,T) (l + J Q (E|6U,nX,)| 4p +E||a(s, n A, 

< C(x 0 ,p,T,K Lip ) ll + j\\\ n X s \\%dsj 



We deduce from Gronwall’s inequality that E|| n A' T ||^ is bounded by a constant which de¬ 
pends only on ccq, P, T and Kl^. 
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Combining this result with m and ED, we get: 


Elinor - nX T fZ < C(xo,p,T,K Lip ) 


1 + log n 


Finally ED becomes: 

n\XT-nX T \\ 2 £ 


< 


C(p) (e||AY - n X T \\% + n\uX T - nXrllS) 


< C(x 0 ,p,T,K Lip ) ( ( 1 + _ bg " V+ [ T E||A S — n X s \\^ds 


And we conclude by Gronwall’s inequality. 
Corollary 3.1. Under Assumption^ 

1, 


Vae[0,-), n a \\X T - n X T \\ z 


0, P—a.s. 


□ 


Proof. Let a £ [0, -). For a p large enough, by Proposition l3T2l we have: 


E 


J2n 2pa \\X T - nX 7 


I 2 V 


7>1 


< 00 


Thus 


and 


^n 2 ^||X T -„A T ||^<oo, P—a.s. 


^>l 


n a \\X T - n X T \\ 


oo —0, P—a.s. 

n—> oo 


□ 


4 Smooth functional approximations for martingales 

Let g : D([0,T],R d ) —> I be a functional which satisfies: 

Assumption 2. g : (D([0, T], R d ), || • Hoc) —> R is continuous with polynomial growth: 
3q € N,3C > 0,Vw £ H([0,T],R d ), \g(cj)\ < C (1 + M *,). 

The (square-integrable) martingale: 

Y(t) = E[g(X T )\T t ] 
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may be represented as a non-anticipative functional of W: 


Y(t) = F(t, Wt) 


where the functional F is square-integrable but fail to be smooth a priori (see [8( for results on 
patlrwise differentiability of F). By Proposition [2751 we have: 


g(X T ) = Y(T) = Y(t) + j X w Y(s) ■ dW(s) 


P—a.s. 


where \7\yY is the weak vertical derivative of Y with respect to W, which cannot be computed 
directly as a pathwise directional derivative unless F is a smooth functional (for example £ 


C&A At))- 


The main idea is to approximate the martingale Y by a sequence of smooth martingales n Y 
which admit a functional representation „Y(s) = F„(s, W s ) with F n £ C^ c ^(At). Then by the 
functional Ito formula, we have: 


/ XMs, W s ) ■ dW{s) = n Y(T) - n Y{t) —* Y(T) - Y(t) = [ V w Y(s) ■ dW{s). 
Jt n ^°° Jt 


One can then use the following estimator for Vifh: 

Z„(s) = V w F„(s,M / s ), 


where the vertical derivative X u F n (s,W s ) = (diF n (s,W s ), 1 < i < d) may be computed as a 
pathwise derivative 


d^s.Ws) 


lim 

h->-0 


Fn(s, Ws + heil[s,T]) 

h 


F n {s, W s ) 


yielding a concrete procedure for computing the estimator. 


We will show in this section that the familiar weak Euler approximation provides a systematic 
way of constructing such smooth functional approximations in the sense of Definition [7| 

Define the concatenation of two cadlag paths w,u/ £ D([0,T],R d ) at time s £ [0, T], which 
we note ui ® u/, as the following cadlag path on [0, T]: 


c o © (jj r — 0J S ® lo' — 

S S 


ui(u) 

U](s) + U}'(u) 


u £ [0, s) 
w'(s) u £ [s, T] 


Observe that: 

V z £ ujg ® u> f = (ca s © u/) © z1.\s,t]- 

Definition 9 (Weak Euler approximation). We define the (level-n) weak Euler approximation 
of F as the functional F n defined by 


F n (s,u] s ) = E 


g( n X(ujs © W T )) 

S 


(13) 
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Applying this functional to the path of the Wiener process W, we obtain a (J r t ) t >o-adapted 
process: 

n Y(a) =F n (s,W s ). 

Using independence of increments of W, we have 


„y(s)=E[ ff ( n X(W T ))|J' s ]=E 


g( n X(W s ®W T ))\T a 

s 


= E 


g{ n X(W s (BB T ))\F S 

S 


where B is any Wiener process independent from W. In particular n Y is a square-integrable 
martingale, so is weakly differentiable in the sense of [7, Theorem 5.8]. We will now show that 
F n is in fact a smooth functional in the sense of Definition [7] 

Theorem 4.1. Under Assumptions |7] and [H the functional F n defined in o is horizontally 
differentiable and infinitely vertically differentiable. 


Proof. First, note that under Assumption [TJ n X(w) is bounded by a polynomial in the variables 
u>(t i—) — w(0), ui(t 2 —) — uj(ti —), • • • , ui(t n —) — ui(t n - 1 —)). Combined with Assumption [2j this 
implies that all expectations in the proof of this theorem are all well defined. 


Let (s,w) G A t with tk < s < tk+i for some 0 < k < n — 1. We start with the vertical 
differentiability of F n at (s,w), which is equivalent to the differentiability at 0 of the map v : 
R d —> R defined by: 


v{z) = F„(s,wf) = E 


g(nX(0J Z s ® Bt)) 

S 


The main idea of the proof is to absorb the dependence with respect to z in the Gaussian density 
function when taking the expectation, which then implies smoothness. 


As we have already shown, „ X(u) z ® Bt) depends only on (u>* ® Br)(ti—) — (w* ® Bt)( 0), 

s s s 

■■■, (wj ® BT){t n —) — (uj z ® BT)(t n -1 —), which are all explicit using the definition of the 

S S 

concatenation. For j < k , we have: 

k 2 © B T ){t j+ 1-) - (w* ® B T )(tj -) = u(t j+ 1-) - 

s s 

In the case where j = fc, we have: 

k © B T ){tk +1 —) — (w, © B T ){tk—) 

S S 

= B(t k +i) - B(s) + uj(s) + z - u(t k -) 

= B(t k+ i) - B(s) + z + u(s) - w(tfe-) 

And for j > k, we have: 

( w s © B T )(t j+1 ~) - (w* ffi B T )(tj—) 

S S 

= B(tj. |_i) - B(s) + w(s) + z— ( B(tj ) - B(s) + oj(s) + z) 

= B(t j+l ) - B(tj) 

Thus we have: 

nX(u z ffi B t ) = p n (oj(ti-) - w(0), • • • ,u(tk~) - u(t k -i-),B(tk+i) - B(s) + z + w(s) - u(tk~), 

s V 

B(tk+ 2 ) — B(tk+ 1 ), • • ■ , B(t n ) — B(t n - 1 )^ 
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where p n : -A/4,n(®0 —> £>([0, T], R d ) is the map defined by ©■ 


Observe from the previous equation that, for a fixed z, the value of n X(tk+i,ui* © Bt) as 

S 

a random variable depends only on a finite number of Gaussian variables: R(4+i) — £?(s), 
B(tk+ 2 ) — -8(4+ 1 ), • • •, B(tj) — B(tj- 1 ). Since the joint distribution of these Gaussian variables 


is explicit, v(z) = E 
dimension. 


© Bt)) 


can be computed explicitly as an integral in finite 


Denoting y = (yi,•• ■ , y n -k) £ A+n-fcQR), 


v(z) = E 
= E 


g( n X T (ujg © Bt)) 

s 

g(p n {u(h-) - +0), • • • ,w(tfc-) - w(tfc-i-), 

8(4+ 1 ) - B(s) + 2 + w(s) - w(tfc-), 8(4+ 2 ) - B{t k+ 1 ), • • • , 8(4) - B{t n - 1 ) 

/ ff(Pn(w(*i-) - w(o),--- ,w(4-) - u(t k -i-),yi + z + w(s) -w(tfe-), 

J(n — k) \ 

n—k 

2/2, - , 2 /n-fc))$(yi,*k+i -s) ^(yi,S)dy 1 dy 2 ---dy n -k 

l =2 

/ ff(p n (w(*l-) - w(0), • • • ,u(tk~) - u{t k -l-),yi + w(s) - w(tfe-), 

J(n — k) \ 

n—k 

2 / 2 , • * * , 2 /n-fc))^(yi -2,4+1 ~ 5 ) II $^)dy\dy2 • • • dyn—k ( 14 ) 


1=2 


with 


$(x,t) = ( 27 r£) 2 exp ( — 


22 


x £ 


the density function of a d-dimensional N(0,tld) variable. Since the only term which depends 
on z in the integrand of is $( 2/1 — z,4 +1 — s), which is a smooth function of z, thus v is 
differentiable at all z G in particular at 0. Hence F n is vertically differentiable at (s,ui) £ A t 
with: for 1 < i < d, 

diF n (s,u) = / g\Pn{u(ti-) - w(0),-■ ■ ,w(4-) - w(tk-i-),yi + w(s) -u(t k -), 

jRdx(n-fe) V 

n—k 

2 / 2 , - - - ,y n -k)) , yi — $(yi,4+i -s) TT $(yi,5)dyidy 2 ---dy n - k 
> 4+i - « 

(-B(4+i) - B(s)) ■ e. 


= E 


g(„X(u) s © Bt))- 


4+i — s 


(15) 


Remark that when s tends towards 4+i, V^i+s, w) may tend to infinity because of the term 
4+i ^ s in the denominator. However in the interval [4,4+i)> V^F n (s,w) behaves well and is 
locally bounded. 

Iterating this procedure, one can show that F n is vertically differentiable to any order. For 
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example, we have: for z G R , 


diF n (s,w z s ) = / g(p n (u(ti-)-u(0),---,u(t k -)-u{tk-i-),yi + z+uj(s)-u(t k -), 

jRdx(n-k) \ 

n—k 

2 / 2 , •• • ,2/n-fe)) , Vl ^ $(2/i,4+i - s) $( 2 //, 5)dyidy 2 • • • cfa/n-fc 

fc+1 s ;=2 


Thus we have: 


d?F n (s,u;) = / g(p„(w(ti-) - w(0), • • • ,w(4-) - w(4-i-),2/i + w(s) - w(4-),2/2,- • • ,2/n-; 

Jjjdx(n-fe) V 

( 77 ^—^ - 7 — -) $( 2 / 1 , 4+1 -s) T| ${yi,S)dy 1 dy 2 ---dy n -k 

V (4+1 - sr 4+1 - s/ ^ 


( / \2 -t \ n—K 

,. yi e% , 2 - 7 -) $( 2 / 1 ,4+i - s ) If $(yi,S)dy 1 dy 2 ---dy n -k 

(4+i - s) 4+i - sy ^ 

And for i ^ j: 

dijF n (s, 10 ) = / g(p„(a;(fi-) - w(0),-• • ,w(4-) - w(4-i-),2/i + w(s) - w(4-), 

jRdx(n-fe) V 

2 / 2 , * * * , 2/n-fc)) ^ fc ^^ /1 g) 2 J ' > ^(2/!, 4+ 1 - s) JJ ${yi,6)dy 1 ---dy n - k 

The horizontal differentiability of F n can be proved similarly. Consider the following map: 
w(h) = F n (s + h, ui s ) = E g( n X(to s ® B T )) , h > 0 

s-\-h 

The objective is to show that w is right-differentiable at 0. 

We assume again that 4 < s < 4 +1 for some 0 < k < n — 1, and we take an h > 0 small 
enough such that s + h < 4 + 1 - Using the same argument and the fact that u> s (s + h) — w(s), 
we have: 

n X(u) a © B t ) = p n {u(ti~) - w(0),- • • ,w(4-) - w(4_i-),B(4+i) - 8(s + ft) + w(s) - w(4-), 

S + ft. 

- 8 ( 4 + 2 ) — 8 ( 4 + 1 ), • • • , 8(4) — 8(4-i)) 

Let y = (yi, ■ ■ ■ ,y n -k) G A4,n-fc(Bi)- We calculate explicitly w(/i): 


tu(/i) = E g( n X T (uj s © 4 t )) 


= E 


.g(p„(w(4-) - w(0),- • • ,w(4-) - w(4_i-),8(4 + i) - 8(s + /i) + w(s) - w(4~), 


8 ( 4 + 2 ) — 8 ( 4 + 1 ), • • • , 8(4) — 8(4-1 


/ g(p n (u;(4-)-w(0),--- ,w(4-)-w(4_i-),j/i + w(s)-w(4-), 

,/R<*X(n-k) \ 

n—k 

2 / 2 ,-•• ,2/n))$(2/i,4+i -s-h) P ${yi,S)dyidy 2 - ■ ■ dy n - k 
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Again the only term which depends on h in the integrand of (1161) is $( 2 / 1 , 4+i — s — h), which 
is a smooth function of h. Therefore F n is horizontally differentiable with: 


VF n (s,ui s ) 





w(0), • • • ,u{t k ~) 



d 

2(4+i 


s) 


\yi \ 2 

2(4 + i - s ) 2 


- u{t k -i-),yi + w(s) -w( 4 -), 

n—k 

®{yi,tk+i ~ s) JJ ^(yi,S)dy 1 ---dy n -k. 
1=2 


□ 


The following result shows that the functional derivatives of F n satisfy the necessary regularity 
conditions for applying the functional Ito formula to F n : 

Theorem 4.2. Under Assumptions Q] and @ F n £ Cif c JA T ). 

Proof. We have already shown in Theorem 14.11 that F n is horizontally differentiable and twice 
vertically differentiable. Using the expressions of VF n , V w F n and Vf,F n obtained in the proof 
of 14.11 and the assumption that g has at most polynomial growth with respect to || • ||oo, we 
observe that in each interval [4,4+ 1 ) with 0 < k < n — 1, VF ni Vand satisfy the 

boundedness-preserving property. We now prove that F n is left-continuous, V u F n and V^.F n 
are right-continuous, and T>F n is continuous at fixed times. 

Let s £ [4,4+i) for some 0 < k < n — 1 and w £ D([0, T], K d ). We first prove that F n is 
right-continuous at (s,ui), and is jointly continuous at (s,w) for s £ (4,4+i)- By definition of 
joint-continuity (or right-continuous), we want to show that: Ve > 0,3 g > 0, V(s',u/) £ A t (for 
the right-continuity, we assume in addition that s' > s), 

d((s,w), (s', w')) <??)=> | F n (s,w) - F n (s',u)')\ < e 


Let (s', oj') £ A t (with s' > s for the right-continuity). We assume that d((s,u), (s', <*/)) < g 
with an g small enough such that s' £ [ 4 , 4 +i) (this is always possible as if s = 4 , we are only 
interested in the right-continuity, thus s' > s). It suffices to prove that | F n (s,u>) — F n (s',u')\ < 
C(s,ui 8 ,g) with C(s,cv s ,i 7 ) a quantity depending only on s, w s and g, and C(s,co s ,g) —> 0. 

77 —>-0 

We use the expression of F n obtained in the proof of Theorem l4.ll Denote y = (y\ , • ■ ■ , y n -k ) £ 
Af+ n _fc(R). We have 


F r 


,(s,w) = / g[Pn(u(tl-) - w(0),-■ • ,w(4-) - w(4-l-),2/l +w(s) - w(4-), 

7Rdx(n-fc) V 

n—k 

2 / 2 ,-•• ,2M))$(2/i,4+i ^ s ) II ${yi 1 S)dyidy 2 ---dy n - k 


and 


F n (s', u/) = [ g(pn{u'(h-) -w'(0),--- ,w'(4-) - w'(4-i-),2/i +w'(s')-u/(4-), 

J^dx(n-fc) V 

n—k 

2 / 2 , * * * , 2 /n))^( 2 / 1 ,4+i s') <&{yi,8)dy\dy 2 - ■ ■ dy n - k . 
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To simplify the notations, we set: 

p(w, s, y ) = p n (v(ti-) - w( 0), • • • , u(t k ~) - w(tfe_i-), yi + w(s) - u(t k -),y 2 , • • ■ , 2/n) 

and 

p(w', s', y ) = p n (w'(ti-) - u/(0), • ■ ■ , u'(t k ~) - w'(4_i-), yi + w'(s') - w'(ifc-), t/ 2 , ■ ■ • , Vn) 
Similarly pt(-) will be the path of p(-) stopped at time i. 

As ||w s — oj s ' ||oo < r), by Lemma 13.11 we have: 

II P(v,s,y) -p{w',s',y) ||oo < C(u> s , y, K Lip , T)rj 

We now control the difference between F n (s,w) and F n (s',u}'). 

\F n {s,uj) - F n {s\u ')| 

« n—/c 

< / |ff(p(w, S, y))$(yi, 4+1 - s) - y(p(w', s', y))$(yi, t k+1 - s')| TT $(yi,6)dy 

J R dx(. n -k ) ^ 

^ (|5(p(w, s, y)) - y(y(w', s', y))|$(yi, t fc+ i - s') 

n—k 

+\g(p{u,s,y))\ • I$(yi,4+1 ^ s) - $(yi,ifc+i - s')Q x $(yi,S)dy 1 dy 2 ■ ■ ■ dy n - k ( 17) 

1=2 



Observe that |$(yi,t fc+ i - s) - $(yi,f fe+ i - s')| < |s - s'| • p{yi,g) < p(yi,y) ■ rj 2 with 


p{ywn) = sup |<9 t $(?/i,f)| 

tG[tk+i-s-r] 2 ,8] 


and we have: 

p(yi,v) 

So the second part of ED can be controlled by: 


—> sup \d t $(yi,t)\= sup 

T?_>0 te[tfc+i-s,<5] te[tk+i-a,s\ 




\yi\ 2 _ d_ 

2 1 2 2 1 


< OO 


,, n—k 

/ \g{p{u,s,y))\-\$(yi,th+i-s)-$(yi,t k +i-s')\ TT $(yi 1 S)dy 1 dy 2 ■ ■ ■ dy n - k < C{s,w 3 ,rj) 

y R t!x(n-fe) J=2 


with 


C(s,u a ,rj) —> 0 . 

77—>0 


For the first part of ED, we use the continuity of g. As y is continuous at p(ui s ,y), we have: 

Iff(p(w, s, y)) - g(p(oj', s', y))| < C(s, y, 77 ) 

with 

C(s,w a ,y,T )) —> 0 

77 —>0 
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and 

$(yi,tfc + i - s') < sup $(y 1 ,t)<oo. 

te[tk+l-s-ri 2 ,S] 

Thus the first part of o can also be bounded by C(s, w s , 77 ) with 

C(s,bJ s ,rj) —> 0. 

ry—>-0 

We conclude that |F n (s,w s ) — F n (s', oj' s ,)\ < C(s,ui s ,r)) depending only on s, w s and 77 , and 
C(s,uj s ,t]) —> 0, which proves the right-continuity of F n and the joint-continuity of F n at all 

77—>-0 

(s,u) £ A t for s / t k , 0 < k < n — 1 . 

The right-continuity of V w F n , F n and the continuity at fixed times of T>F n can be similarly 
shown from the expressions of V \jF n , V^,F„ and DF n obtained in Theorem ld.il Now it remains to 
show that for 1 < k < n and ui £ £>([0, T], R d ), F n is left-continuous at (t k ,w)- Let (s', <*/) £ A t 
with s' < t k such that d((t k ,uj), ( s',w')) < 17 . We choose an 77 small enough in order that s' £ 
[tk-i,tk), and we want to show that \F n (t k ,u) - F n (s',w')\ < C(t k , u>t k ,v) f° r some C{tk,ut k ,i l) 
depending only on t k , uit k and 77 with C(t k ,uit k ,r]) —► 0- 

77^0 

We first decompose \F n (t k ,uj) — F n (s' ,uj')\ into two terms: 

| F n (t k ,u) - F n (s',u ;')| < | F n (t k ,u) - F n (s r , ov)| + \F n (s',uj s >) - F n (s',uj')\ 

For the second part, as F n is continuous at fixed time s' by the first part of the proof, and 

||t*V - w'/||oo < Vi we have |F„(s',uv) -F„(s',w')| < C(t k , u tk , 77 ) with C(t k , u} tk , 7?) —> 0. 

77—^0 


For the first part | F n (t k ,uj) — F n (s', uv)|, the difficulty is that s' and t k no longer lie in 
the same interval, thus we need to perform one more integration for F n (s',oj s i) compared to 
F n (t k , u). Let y = ( 7 / 1 , • • • , y„- k ) £ M d ,n-k( and y' £ R d . Using again the expression of F n 
we have obtained in the proof of Theorem 14. II we have: 


F n (t k ,w) — 


/ g\P n {u(ti-) -w(0),--- ,u}{t k -)-u(th-i-),yi+u{th)-u(tk-), 

jR^x(n-fc) V 

n—k 

2 / 2 ,-•• ,2/n-k)) dyi®{yi,6) 


and 


F n (s',uv) = / ff(Pn(w(ti-) - w(0),--• ,y' +uj(s')-uj{t k - 1 -),y 1 ,--- ,y n )\ 

J^dx(n-k+l) \ / 

n—k 

®(y', tk - s')dy' dyi$(yi,6) 

1=1 

= / (/ ff(p(w(ti-)-w(0),--- ,y' +w(s')-u(t k - 1 -),yi,---,y n )) 

JR d V jRdx(n-fc) 

n—k 

n ®{yi,8)dyt^${y',t k - s')dy' 


7=1 
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We now define £ : R d —> K by: for y' £ K. d , 


C(y') = f g(Pn(u(ti~) - w(0), • • ■ , y' + w(s') 

J ]RdX(n-fc) 


n—k 

w(ffe-i-),yi, • • • >2/n)) ${yi,S)dyi 

1=1 


By Lemma 13.11 and the continuity of g with respect to || • ||oo, the map 

y' g(p(v{ti~) - w(0), • • • , y' + w(s') - w(ffe-i-), yi, ■ ■ ■ , y n )) 

is continuous. As g has at most polynomial growth with respect to || • Hoo, by the dominated 
convergence theorem, £ is also continuous. Moreover, £ has at most polynomial growth. And as 
ffc — s' < rj 2 , we have 


F n (s',ui s >) 


/ C(2 /)<S>(y',t k -s')dy' 

JR d 

[ (ay')-mMy\tk-s'w+m- 

JR d 


with 

/ (Cfo')-C(o)W,tk 

JRd 

and C(tk,u>t k ,V) —> 0. 

77—>■ 0 



< C(t k ,u tk ,r]) 


It remains to control the difference between F n (tk,u) and £(0). We remark that: 

„ n—k 

C(0) = / g(p n (u{ti-) -w(0),--- ,u(s') -w(4_i-),yi,--- ,y n )) TT ®(yi,6)dyt 

jRdX(n-k) ; = 1 

= Fn(J'k^ C^s') 


= E 


S(h^t( w s' © -Bt)) 

tfc 


As ||uv - wtJloo < ||uv - w',||oo + ||w tfc - w',||oo < 2 1 ), again by the continuity of F n at fixed 
time tk established in the first part of the proof, we have: 


with C(tk,u>t k , rf) —> 0. 

77—>-0 

We conclude that 


I F n {t k ,uj) - C(0)| < C(tk,u>t k , rf) 


I F n (t k ,u) - F n (s',uj')\ < C(tk,u} tk ,r)) 


with C(tk,u>t k ,ri) —t 0, which proves the left-continuity of F n at (tk,u>). 


77—>•() 


Corollary 4.1. Under Assumptions [7] for any t G [0,T) we have: 

F n (T,W T ) - F n (t,W t ) = J V u F n (a,W.) ■ dW(s), P - a.s. 


□ 


(18) 


Proof. As F n £ ^locri At), we can apply the functional Ito formula Proposition [2721 and we 
remark that the finite variation term is zero as n Y(s ) = F n (s , W s ) is a martingale. □ 
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Remark 4.1. We can also verify using directly the expressions we have obtained in Theorem \4-l\ 
for T>F n and X^F n that the finite variation terms in HID cancel each other. By the functional 

ltd formula, the finite variation term in CEE]) equals to VF n (s,W s ) + -tr(VlF n (s,W s )). And for 
(s,w) £ At with s £ [4,4+ 1 ), we have: 


te{VlF n (s,u)) ='}Td?F n (s,uj) 


I 


R dx(n-k) 


g(p n (uj(t 1 -) - w(0), • • • ,w(tfc-) - + w(s) - U)(tk-),y2,- ■ ■ ,Vn)) 


a / ( \2 ^ \ Ti—K 

( (+ yi e ' \2 - 7 - J - s) n $( yi,5)dy 1 • •• dy n - k 

\ (tk +1 - s) 2 4+i - sj fj* 2 

/ g(Pn(u(t 1 -) - w(0), • • • ,w(tfe-) - w(4_i-), J/1 + w(s) - uj{t k -),y 2 ,-■ ■ ,y n )) 

J R<*X(n-fc) 

( \yi\ 2 


d \ Fi—k 

-) $( 2 / 1 ,4+1 - s) TT ^{yu$)dyi ■ ■ ■ dy n -k 

1 - f=2 


V(4+1 - s) 2 4+1 

= -2 VF n (s,uj s ) 

which confirms that F n is a solution of the path-dependent Kolmogorov equation Q Sec. 5]: 

VF n {s , Ws) + ^tr (ylF n (s, W,)) = 0. 


5 Convergence and error analysis 


In this section, we analyze the convergence rate of our approximation method. After having 
constructed a sequence of smooth functionals F n ( Theorem 14.11 and Theorem 14.21) . we can now 
approximate Vi+F by: 

Zn(s) = V w F„(s, W s ) 

which, in contrast to the weak derivative VwY", is computable as a pathwise directional deriva¬ 
tive. In practice, V u F n (s, W s ) can be computed numerically via a finite difference method or a 
Monte-Carlo method using the expression m of v w f„. 


For t £ [0,T], the quantity of interest is the integral of Vn/h — Z n along the path of W 
between t and T, i.e. 


[ (V W Y ~Z n )-dW= [ V w Y(s) ■ dW(s) - [ V„F n (s, W.) ■ dW(s ) 

Jt Jt Jt 

By the martingale representation formula Proposition 12.31 and Corollary 14.11 we have P-a.s. 
J iy W Y — Z n ) ■ dW = Y(T) — Y(t) — ( n Y(T) — n Y(t)) 


= g(X T )~g(nX T (W T ))-E 


g(X T ) - g( n X T (W t ® B T ))\F t 
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where n X is the path of the piecewise constant Euler-Maruyama scheme defined in 0. Remark 
that by definition of the concatenation operation, if B and W are two independent Brownian 
motions, we have: 


E 


g{ n X T (W t ®B T ))\Tt 


= E 


g(nX T (W t ®w T ))\Xt 


= E[g( n X T (W T ))\Xt] 


Proposition 5.1. Under Assumptions^ and\^ 


Vf G [0,T], [ (y w Y-Z n )-dW-^ 0, 

Jt n ^°° 




Proof. We have already shown that: 

(V w r - Z n ) ■ dX = g{X T ) - g{ n X T ) - E [g(X T ) - g{ n X T )\F t ] 

As g is continuous with respect to || • Hoo, by Corollary 13.11 we have: 

g{X T )~g( n X T ) —s> 0, P-a.s. 

n—too 



Under assumption [2] the sequence of martingales Y n = E[g( n X T )\T u is bounded in L 2 therefore 
uniformly integrable, hence the result. □ 

Corollary 5.1. Under Assumptions^ and\^ 

r T 


Vf g [o,T], 


(V W Y - Z n ) ■ dW 


2 V 


Vp > 1 


Proof, g has at most polynomial growth with respect to || • ||oo, which ensures the uniform 
integrability of \g( n X T )\ 2p . We conclude by applying the dominated convergence theorem. □ 


Under a slightly stronger assumption on g we can obtain a rate of convergence: 

Theorem 5.1 (Rate of convergence). Let p > 1 and assume g : (D([0, T], R d ), ||.||oo) —> R is 
Lipschitz- continuous: 


3gLi P >0, Vw,u/ G £>([0,T],R d ), \g(u) - g{u')\ < g Lip ||w - u/||oo- 

Under Assumptions |7J the L 2p -error of the approximation Z n of V\yY along the path of W 
between t and T is bounded by: 


E 



z n ) ■ dW 


2p\ 


^ U(xo , p, T, E-Lip, gLip) 


^1 + logre ^ p 


Vp > 1 


where the constant C depends only on Xq,p,T, Kli P and gLip- In particular: 


Va G [0, -), 



Z n ) ■ dW 


P—a.s. 
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Proof. This result is a consequence of Proposition 13.21 since 


(X w Y(s) - z n {s)) ■ dW (s) 


2 P 


< 

< 

< 


II 9(Xt) - g(nX T )h P + m\g(X T ) - g( n X T )\P t ]\\ 2p 

2\\g(X T ) - g( n X T )\\ 2p 

ZgLipW sup |X(s) - n x(s)||| 2p . 

«e[o,T] 


□ 


The following example how our result may be used to construct explicit approximations with 
controlled convergence rates for conditional expectation of non-smooth functionals: 

Example 5.1. Let 

g{u) = ip(w{T), sup ||w(i)||) 
te[o,T] 

where if € x R + ,K) is a continuous function with polynomial growth, and set Y(t) = 

E[g(XT)\Xt\. Then g satisfies Assumption^ and our approximation method applies. Moreover, 
if ip is Lipschitz-continuous, then Theorem \5.1\ yields an explicit control of the approximation 
error with a ydog n/n bound. 


6 Comparison with approaches based on the Malliavin cal¬ 
culus 


The vertical derivative S7xY(t) which appears in the martingale representation formula may be 
viewed as a ’sensitivity’ of the martingale Y to the initial condition X(t). Thus, our method is 
related to methods previously proposed for ’sensitivity analysis’ of Wiener functionals. 

One can roughly classify such methods into two categories |2|: methods that differentiate 
paths and methods that differentiate densities. When the density of the functional is known, 
the sensitivity of an expectation with respect to some parameter is to differentiate directly the 
density function with respect to the parameter. However, as this is almost never the case in a 
general diffusion model, let alone a non-Markovian model, alternative methods, are used: these 
consist of differentiating either the functional g or the process with respect to the parameter 
under the expectation sign, then estimating the expectation with the Monte-Carlo method. To 
differentiate process, one required the existence of the so-called first variation process, which 
requires the regularity of the coefficients of the SDE satisfied by X. 

Sensitivity estimators for non-smooth functionals may be computed using Malliavin calculus: 
this approach, proposed by Fournie et al. !18| and developed by Cvitanic, Ma and Zhang [5], 
Fournie et al. E|, Gobet and Kohatsu-Higa Kohatsu-Higa and Montero [25], Davis and 
Johansson E and others, uses the Malliavin integration-by-parts formula on Wiener space in 
the case where g is not smooth. These methods require quite demanding regularity assumptions 
(differentiability and ellipticity condition on er for example) on the coefficients of the initial SDE 
satisfied by X. 
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By contrast, the approximation method presented here allows for any continuous functional g 
with polynomial growth and requires only a (Lipschitz-)continuity assumption on the functional 
a and allows for degenrate coefficients. It is thus applicable to a wider range of examples than the 
Malliavin approach, while being arguably simpler from a computational viewpoint. Contrarily 
to the Malliavin approach which involves differentiating in the Malliavin sense, then discretizing 
the tangent process, our method involves discretizing then differentiating (the Euler scheme) 
which, as argued in J2], has its computational advantages. 

In our setting, we have F n £ Cj 1 ^ (A t) which is sufficient for obtaining an approximation of 
martingale representations via the functional Ito formula. One can ask if the Euler approximation 
n X can also be used to obtain a Clark-Haussmann-Ocone type formula, and in this case, whether 
the pathwise vertical derivative V w .F n (i, Wt) leads to the same representation as the Clark- 
Haussmann-Ocone formula. 

For n £ N, define H n = g( n X T (WT )) with n X the weak piecewise constant Euler-Maruyama 
scheme defined by (ED. By the definition of n X, the random variable H n actually depends only 
on a finite number of Gaussian variables: W(tx), Wfo) — W(tx), ■ ■ ■, W(t n ) — W(t n -i), thus it 
can be written as: 


H n = h n (W(tx), W(t 2 ) - • • , W(t n ) - W{t n _x)) 

with h n : Ald, n (R) —> M (h n is actually g op with p defined by ©)• 

Clearly if h n is a smooth function with polynomial growth, then H n £ D 1 ' 2 with Malliavin 
derivative EZj: 

H>tH n = ( H>tH n , 1 < k < d) £ R d , where 

n—1 

B $H n = Y l d kj h n {W{t 1 ),W{t 2 )-W{t 1 ),--- ,W(t n ) t£[ 0 ,T) 

3=0 

In this case, assume that t £ [tj,tj+ 1 ) for some 0 < j < n — 1, we have: for 1 < k < d, 
E[B*H n \T t \ 

= E [d kj h n (W(tx), W{t 2 ) - W{tx), ■ • • , W{t n ) - W(t n -x))\T t ] 

= E[dkjh n (u(tx),u(t2) -u(tx),--- ,u(tj) -w(t j - 1 ),W(t j+ 1) - W(t)+u}(t) 

W(t j+2 ) - W(t j+1 ), • • • , W(t n ) - W(t n -x))] L=W t 

Jd_ 

idh' 

W(t j+2 ) - W(t j+ 1 ), • • • , W(t„) - W(t n - i))U = 0 
= ^(E[h n (u](tx),u(t2) - u(tx),-- ■ - uj(tj-x),W{t j+1 ) - W(t)+ui(t) - u(tj) + he k , 

W(tj +2 ) - W(t j+1 ), ■■■, W(t n ) - W(t n .x))] Ut=Wt ) \h=0 

which is none other than 

Fn{t, Wt + he k l[t,T]) — Fn{t, Wt) 


= E ^/i n (w(ti),w(t 2 ) - w(H), • • • - oj(tj-x),W(t j+ x) - W{t) +w(f) - u(tj) + he k , 


I u t =Wt 


d k F n (t, W t ) = lim 

h—> 0 


h 


So in the case where h n are smooth, our method provides the same result as given by the Clark- 
Haussmann-Ocone formula applied to h n . However, in our framework, as the functional g is only 
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assumed to be continuous with polynomial growth, the function h n may fail to be differentiable. 
So, even in the cylindrical case, it is not clear whether the random variable H n is differentiable 
in the Malliavin sense, and even if it is the case, it is difficult to obtain an explicit form for 
E using the Malliavin calculus. 

However, our approximation method applies even in the cases where H n is not differentiable 
in the Malliavin sense: indeed, as shown in Section 3, as soon as H is square-integrable, the 
martingale n Y(t) = E[H n \Tt] has a smooth functional representation which is differentiable 
in the pathwise sense, even though H n is not differentiable, neither in a patlrwise nor in the 
Malliavin sense. 


References 

[1] Jean-Michel Bismut. A generalized formula of Ito and some other properties of stochastic 
flows. Z. Wahrsch. Verw. Gebiete , 55(3):331-350, 1981. 

[2] Nan Chen and Paul Glasserman. Malliavin Greeks without Malliavin calculus. Stochastic 
Processes and their Applications, 117( 11): 1689—1723, 2007. 

[3] J. M. C. Clark. The representation of functionals of Brownian motion by stochastic integrals. 
Ann. Math. Statist., 41:1282 1295, 1970. 

[4] Rama Cont. Functional Ito calculus and path-dependent Kolmogorov equations (Lecture 
Notes of the Barcelona Summer School on Stochastic Analysis, July 2012). Advanced 
Courses in Mathematics - CRM Barcelona. 2012. 

[5] Rama Cont and David-Antoine Fournie. Change of variable formulas for non-anticipative 
functionals on path space. J. Fund. Anal, 259(4):1043-1072, 2010. 

[6] Rama Cont and David-Antoine Fournie. A functional extension of the Ito formula. C. R. 
Math. Acad. Sci. Paris, 348(l-2):57-61, 2010. 

[7] Rama Cont and David-Antoine Fournie. Functional Ito calculus and stochastic integral 
representation of martingales. The Annals of Probability, 41(1):109-133, 2013. 

[8] Rama Cont and Candia Riga. Pathwise analysis and robustness of hedging strategies for 
path-dependent derivatives. Working paper, 2014. 

[9] Jaksa Cvitanic, Jin Ma, and Jianfeng Zhang. Efficient computation of hedging portfolios 
for options with discontinuous payoffs. Mathematical Finance, 13(1):135—151, 2003. 

[10] M. H. A. Davis. Functionals of diffusion processes as stochastic integrals. Math. Proc. 
Cambridge Philos. Soc ., 87(1):157 166, 1980. 

[11] Mark H. A. Davis and Martin P. Johansson. Malliavin Monte Carlo Greeks for jump diffu¬ 
sions. Stochastic Process. Appl, 116(1):101129, 2006. 

[12] Bruno Dupire. Functional Ito calculus. Bloomberg Portfolio Research paper, (2009-04), 2009. 

[13] Robert J. Elliott and Michael Kohlmann. A short proof of a martingale representation 
result. Statist. Probab. Lett., 6(5):327-329, 1988. 


26 



[14] Olivier Faure. Simulation du mouvement brownien et des diffusions. PhD thesis, Ecole 
Nationale des Ponts et Chaussees, 1992. 

[15] P. J. Fitzsimmons and B. Rajeev. A new approach to the martingale representation theorem. 
Stochastics , 81(5):467-476, 2009. 

[16] David-Antoine Fournie. Functional Ito calculus and applications. ProQuest LLC, Ann Arbor, 
MI, 2010. Ph.D. Thesis-Columbia University. 

[17] Eric Fournie, Jean-Micliel Lasry, Jerome Lebuchoux, and Pierre-Louis Lions. Applications 
of Malliavin calculus to Monte-Carlo methods in finance. II. Finance Stock., 5(2):201 236, 
2001 . 

[18] Eric Fournie, Jean-Michel Lasry, Jerome Lebuchoux, Pierre-Louis Lions, and Nizar Touzi. 
Applications of Malliavin calculus to Monte Carlo methods in finance. Finance and Stochas¬ 
tics, 3(4):391-412, 1999. 

[19] Emmanuel Gobet and Arturo Kohatsu-Higa. Computation of Greeks for barrier and look- 
back options using Malliavin calculus. Electron. Comm. Probab., 8:51-62 (electronic), 2003. 

[20] Emmanuel Gobet and Stephane Menozzi. Discrete sampling of functionals of Ito processes. 
In Seminaire de Probabilites XL, volume 1899 of Lecture Notes in Math., pages 355-374. 
Springer, Berlin, 2007. 

[21] U. G. Haussmann. On the integral representation of functionals of Ito processes. Stochastics, 
3(1):17-27, 1979. 

[22] Desmond J. Higharn, Xuerong Mao, and Andrew M. Stuart. Strong convergence of Euler- 
type methods for nonlinear stochastic differential equations. SIAM J. Numer. Anal., 
40(3): 1041-1063, 2002. 

[23] Jean Jacod, Sylvie Meleard, and Philip Protter. Explicit form and robustness of martingale 
representations. Ann. Probab., 28(4):1747 1780, 2000. 

[24] Ioannis Karatzas, Daniel L. Ocone, and Jinlu Li. An extension of Clark’s formula. Stochastics 
Stochastics Rep., 37(3):127 131, 1991. 

[25] Arturo Kohatsu-Higa and Miquel Montero. Malliavin calculus in finance. In Handbook of 
computational and numerical methods in finance, pages 111-174. Birklrauser, 2004. 

[26] Dorival Leao and Alberto Ohashi. Weak approximations for Wiener functionals. Ann. Appl. 
Probab., 23(4):1660-1691, 2013. 

[27] David Nualart. Malliavin calculus and its applications, volume 110 of CBMS Regional 
Conference Series in Mathematics. American Mathematical Society, Providence, RI, 2009. 

[28] Daniel Ocone. Malliavin’s calculus and stochastic integral representations of functionals of 
diffusion processes. Stochastics, 12(3-4):161 185, 1984. 

[29] Gilles Pages. Introduction to numerical probability for finance. http: 
//www.proba.jussieu.fr/pageperso/pages, 2014. 

[30] E. Pardoux and S. Peng. Backward stochastic differential equations and quasilinear parabolic 
partial differential equations. In Stochastic partial differential equations and their applica¬ 
tions (Charlotte, NC, 1991), volume 176 of Lecture Notes in Control and Inform. Sci., pages 
200-217. Springer, Berlin, 1992. 


27 



[31] E. Pardoux and D. Talay. Discretization and simulation of stochastic differential equations. 
Acta Applicandae Mathematica, 3(1):23—4T, 1985. 

[32] Philip E. Protter. Stochastic integration and differential equations, volume 21 of Stochastic 
Modelling and Applied Probability. Springer-Verlag, Berlin, 2005. Second edition. 

[33] Daniel Revuz and Marc Yor. Continuous martingales and Brownian motion , volume 293 
of Grundlehren der Mathematischen Wissenschaften. Springer-Verlag, Berlin, third edition, 
1999. 


28 



